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Abstract: We present a convenient method to generate vector beams 
of light having polarization singularities on their axis, via partial spin-to- 
orbital angular momentum conversion in a suitably patterned liquid crystal 
cell. The resulting polarization patterns exhibit a C-point on the beam axis 
and an L-line loop around it, and may have different geometrical structures 
such as "lemon", "star", and "spiral". Our generation method allows us to 
control the radius of L-line loop around the central C-point. Moreover, we 
investigate the free-air propagation of these fields across a Rayleigh range. 
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1. Introduction 

It has been widely studied and directly observed that both longitudinal and transverse waves 
may have structures rich of singular points [fl]|2]. Optical light beams with singular scalar and 
vectorial features received a particular attention [2"j. In the scalar case, these singularities in 
the beam transverse plane are known as "wave dislocations" or "optical vortices" 0][3]. They 
correspond to points where the optical field amplitude vanishes and the optical phase becomes 
undefined. Near the singular point, the optical phase changes continuously with a specific topo- 
logical charge, given by the phase change in units of 2n along a closed path surrounding the 
singularity 0. Beams with optical vortices carry orbital angular momentum (OAM), which is 
simply proportional to the topological charge when the vortex is located on the beam axis J4). 

When also the vectorial nature of the optical field is taken into account, the variety of pos- 
sible singular points in the transverse plane is increased Q UID . In the paraxial approxima- 
tion, the optical field is purely transverse, i.e. orthogonal to the propagation direction. In each 
point of the transverse plane, the direction of the electric field changes in time tracing an el- 
lipse, whose ellipticity and orientation define the local "polarization" state of the field. The 
polarization can be also parametrized by introducing the three reduced Stokes parameters, s\, 
S2 and S3, or the polar and azimuthal angles on the Poincare sphere. The vector-field singu- 
larities form points or lines in the transverse plane where one of the polarization parameters 
is undefined. Examples are the so-called C-points and L-lines, where the orientation and the 
handedness of the polarization ellipse are undefined, respectively 0171. Optical beams exhibit- 
ing vector singularities of this kind have been given various names in the literature, such as 
for example (full)-Poincare beams II121I13I . Here, we simply call them "polarization-singular 
beams'" (PSB). The PSB may have different transverse polarization patterns, characterized by a 
different integer or half-integer topological index rj as defined by the variation of the polariza- 
tion ellipse orientation for a closed path around a C-point, such as "star", "monstar", "lemon", 
"spiral" (with a central "source/sink" singularity), "hyperbolic" (with a central "saddle" singu- 
larity), etc. Some examples are shown in Fig. [I] Optical PSB possessing high values of positive 
and negative topological charge form so-called "polarization flowers" and "hyperbolic webs", 
respectively |8l . 

Polarization singularities are typically generated either by interferometric methods (see, 
e.g., lfT4l ). which is usually very delicate and unstable, or by propagation through suitable 
birefringent solid crystals (see, e.g., lfj"5l and references therein), a method which does not lend 
itself to a flexible control of the output. In this paper, we study the generation and propagation 
dynamics of vectorial polarization singular beams generated by "q-plates" lfl6l . i.e. birefrin- 
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Fig. 1. Polarization patterns of three different polarization singular beams: "star" (index 
T] = —1/2), "lemon" (ij = 1/2) and "spiral" (rj = 1). The central point at the origin in- 
dicates the C-point, while dashed blue lines show the L-line circles with radius po- The 
handednesses of the polarization ellipses in the outer region, p > po, and in the inner one, 
p < Po, are opposite. 



gent liquid crystal cells having a singular pattern of the optical axis. In a recent work, we have 
reported the generation of singular vortex beams with uniform polarization ellipticity and a 
central vortex singularity with vanishing intensity, as generated by "tuned" q-plates fTD . Here, 
by exploiting "untuned" q-plates, we conveniently generate and control beams having a space- 
variant polarization ellipticity, with a central C-point singularity, and a lemon, star or spiral 
ellipse-orientation pattern. Finally, we analyze the varying polarization structures in different 
longitudinal positions along the beam, so as to investigate the effect of propagation. 

2. Polarization singular beams 

A class of polarization singular beams can be generated by the superposition of two optical 
fields bearing orthogonal polarization states and different OAM values. Here, we consider in 
particular a beam described by the following expression: 

|PSB) =cos (^j LG ,o(p^,z)\L)+e ,a S m (|j LG , m (p,(j),z)\R) (1) 

where p,(j),z are cylindrical coordinates, 8 is a tuning parameter, \L) and \R) stand for left- and 
right-circular polarizations, respectively, a is a phase, and LG p . m denotes a Laguerre-Gauss 
mode of radial index p and OAM index m, with a given waist location (z = 0) and radius wq. 

The LG modes with m ^ have a doughnut shape, with vanishing intensity at the center. 
The PSB resulting from Eq. (fl~|i then have a polarization pattern, shown in Fig.Q] characterized 
by the following four elements: (i) at the center (that is the beam axis), the polarization is left- 
circular (C-point), since the last term in Eq. (Q3 vanishes here; (ii) as the radius increases, the 
second term contributes and changes the circular polarizations into left-hand elliptical; (iii) at a 
certain radius p = po(5), where the two waves in Eq. (Q]i have equal amplitude, the polarization 
becomes exactly linear (L-lines); (iv) at larger radii, that is for p > po, the second term in Eq. (Q~|i 
prevails and the polarization state changes into right-handed elliptical. It is worth noting that, 
in this case, the polarization pattern has a topological index rj = m/2. The initial angle at 
azimuthal location <p = of the polarization ellipse major axis is fixed by the relative phase a 
in Eq. (Q~|). Of course, the \L) and \R) polarizations in Eq. fl} can also be swapped, leading to 
an inverted polarization ellipticitiy everywhere. 
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Fig. 2. Experimental setup. A He-Ne laser beam is first spatial-mode filtered by focusing 
via a microscope objective (4X) into a 50 flm pinhole and recollimating by a lens (f). Its 
polarization is then prepared in a left (right) circular state by a polarizing beam-splitter 
(PBS) and a quarter-wave plate (Q). An untuned q-plate then transforms the beam into 
a PSB. The PSB polarization pattern was analyzed by a quarter- wave plate, a half- wave- 
plate (H) and another polarizer, followed by imaging on a CCD camera. The S parameter of 
the PSB was adjusted by electrically tuning the q-plate. In order to study the propagation 
dynamics of the PSB, the CCD camera was mounted on a translation stage and moved 
around the focal plane of an output lens. Upper insets show the intensity and reconstructed 
polarization patterns, in the near field, of the lemon, star, and spiral beams generated by 
q-plates with q = 1 /2, — 1/2, 1, respectively. 

We generate PSB by exploiting partial spin-to-orbital angular momentum conversion in q- 
plates. The optical action of q-plates is described in HMHrSIi . Depending on its birefringent 
retardation 8, the q-plate converts a varying fraction of the spin angular momentum change 
suffered by the light passing through the device into OAM with m = ±2q, where q is the q-plate 
topological charge and the sign is determined by the input circular polarization handedness. The 
conversion is full when the q-plate optical phase retardation is 8 = n (or any odd multiple of k). 
For other values of 8, the q-plate is said to be untuned and the spin-to-orbital conversion is only 
partial. In the untuned q-plate a fraction of the incoming beam passes through unchanged. The 
field at the exit plane of an untuned q-plate for a Gaussian left-circular polarized input beam is 
given by the following expression: 

U-|L) =HyGG_ |9|VW5 (p,d/n) 

where HyGG ^^/j t/v /2 denotes the amplitude profile of a hypergeometric-Gaussian beam lETTl . 
describing with good approximation the output mode profile of a q-plate of charge q, thickness 
d and average refractive index h for a Gaussian input (see Ref. 0221 for more details), and the 
phase a in the second term depends on the initial angle (on the x axis) of the q-plate optical- 
axis pattern 01611181 . The converted and non-converted terms in Eq. (|2}, that is in the near field, 
have a common amplitude profile, which can therefore be factorized. The far-field profiles 
produced by the converted and non-converted parts of the beam are different, however, because 
the different OAM values affect the beam propagation. The HyGG term appearing in Eq. (O, 
after multiplication by the phase factor exp (2iq<j>), can be expanded in a series of LG modes 
iBTl . Truncating this expansion to the first useful order (which is a good approximation in the 
far field), Eq. (01 is reduced to Eq. (HJ with m = 2q. The optical retardation 8 of the q-plate, 
that is the q-plate tuning, can be adjusted continuously by applying an appropriate voltage [20J. 
Different polarization topologies can be obtained by using q-plates with different topological 
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Fig. 3. Intensity distribution and reconstructed polarization patterns of beam generated by a 
q-plate with q — 1/2 for seven different optical retardations: (a) 8 = (or 2k), (b) 8 = 71 / 4, 
(c) 8 = k/2, (d) 8 = 3k/4, (e) 8 = %, (f) 5 = 5k/4, and (g) 8 = 3n/2, The corresponding 
L-line radii relative to the beam waist wo are the following: (a) undefined (b) po = 2.0wq, 
(c) po = l-4wo, (d) po = I.Iwq, (e) undefined, (f) po = l.Owo, and (g) po = 1.5wo- 

charges. In this work, we used q = —1/2, +1/2, and +1, generating correspondingly equal 
indices of the polarization patterns. 

3. Experiment 

The experiment layout is shown in Fig. [2] The generated PSB polarization distribution was stud- 
ied both in the near and far fields generated in the image and focal planes of a lens. The polar- 
ization distribution was analyzed by polarization tomography: the local Stokes parameter maps 
have been calculated from the intensity images taken by a CCD camera (Sony AS-638CL), after 
projection onto horizontal (H), vertical (V), anti-diagonal (A), diagonal (D), L and R polariza- 
tion states, respectively. The polarization pattern was thus entirely reconstructed, as shown in 
form of little ellipses in Fig. |2]and subsequent figures. In order to reduce noise and experimen- 
tal errors, the average intensity was recorded in a grid of 20x20 pixel square area 0171 . The 
missing ellipses in the experimentally reconstructed polarization patterns (e.g., in Fig. [3]l cor- 
respond to regions with low intensity, so that the polarization reconstruction was too affected 
by stray light. The superposition coefficients in Eq. <J2J have been adjusted by changing the 
q-plate optical retardation 8 with the applied voltage. The radius where the L-line forms in the 
polarization pattern depends on the ratio between the absolute value of the two coefficients in 
Eq. (O and, hence, on the voltage applied to the q-plate. In particular, for optical retardation 
8 = and 8 = K the beam profile is Gaussian fully left-circular and doughnut-shaped fully 
right-circular, respectively. In the intermediate case, the pattern is a superposition and a circular 
L-line appears in the polarization pattern, whose radius po depends on 8. From Eq. ([TJ one finds 
Po = w o/ (V2 tan 8/2). Figure © shows the intensity and reconstructed polarization pattern for 
different q-plate retardations in steps At) = n/4, and the corresponding measured L-lines radii. 
As shown in the figure, also the intensity profile changes from a Gaussian (left-circular) to a 
doughnut (right-circular) shape. In Fig.|3](c) and (g), the optical retardations are n/2 and 3n/2, 
respectively; the n difference in the relative phases results into a 180° rotation of the polariza- 
tion patterns. The topological index of the polarization pattern, however, depends only on the 
g-value of the q-plate, and in the case of Fig. [3] we have X\ = q = 1/2. 

In a second experiment, we studied the dynamics of different PSBs with rj = — 1/2, +1 /2 
and +1 topologies under free-air propagation for fixed parameter 8 = n/2. To this purpose, 
we moved the CCD camera along the propagation axis around the imaging-lens focal plane. 
We recorded the beam polarization patterns at six different planes in the range — zr < z < Zr, 
where zr is the lens Rayleigh parameter and z = corresponds to the beam waist location. The 



Fig. 4. Reconstructed experimental polarization patterns of different PSB beams, (a) 
Tj = m/2 = +1/2, (b) 77 = m/2 = — 1/2 and (c) TJ = m/2= +1. Patterns have been recon- 
structed by measuring the maps of reduced Stokes parameters in six different longitudinal 
planes within the beam Rayleigh range, from — zr to +zr- The corresponding rotation of 
polarization patterns for (a) and (b) are 90° and 88° = (30 + 30 + 28)°, respectively. 

experimental results are shown in Fig. |4] As it can be seen, the polarization pattern evolves 
during propagation. Indeed LGoo and LGo.,„ modes of Eq. (fTJ have different z dependences of 
their Gouy phases; the relative Gouy phase between them is given by y/ = \m\ arctan (z/zo), 
which in the explored region varies in the range \m\[— n/A,n/A\. In the case m = ±1 (i.e. 
77 = 1 /2) the phase evolution leads to a rigid rotation of the whole polarization structure by an 
angle equal to \jf [as shown in Fig.|4](a) and (b)], while when m = 2 (i.e. tj = 1), this dephasing 
leads to a more complex pattern dynamics, in which the polarization structure changes from 
radial to spiral and then to azimuthal (Fig.|4](c)). 

4. Conclusion 

We reported a non-interferometric technique to generate beams carrying a polarization singu- 
larity on their axis, based on a partial spin-to-orbital angular momentum conversion in a de- 
tuned q-plate device. By this approach, we generated beams with C-points on their optical axis, 
surrounded by a circular L-line whose radius can be controlled by the voltage applied to the 
q-plate. We demonstrated different patterns such as stars, lemons, and spirals, corresponding to 
vector topological indices of — 1 /2, 1 /2 and +1, respectively. We have also studied the pattern 
evolution in free-space propagation. Such polarization patterns could find practical application 
in lithographic optical techniques, as for the examples discussed in Refs. 112311241 . 
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